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Summary. The detailed theoretical understanding of quantum spin dynamics in various 
molecular magnets is an important step on the roadway to technological applications of these 
systems. Quantum effects in both ferromagnetic and antiferromagnetic molecular clusters 
are, by now, theoretically well understood. Ferromagnetic molecular clusters allow one to 
study the interplay of incoherent quantum tunneling and thermally activated transitions 
between states with different spin orientation. The Berry phase oscillations found in Fes 
are signatures of the quantum mechanical interference of different tunneling paths. Anti- 
ferromagnetic molecular clusters are promising candidates for the observation of coherent 
quantum tunneling on the mesoscopic scale. Although challenging, applications of molecu- 
lar magnetic clusters for data storage and quantum data processing are within experimental 
reach already with present day technology. 

Introduction 

Molecular magnets have attracted considerable interest recently because of their potential 
for data storage and data processing 0. In addition to possible future technological ap- 
plications, molecular magnets are also interesting from an academic point of view because 
they show quantum effects on the mesoscopic scale || in the form of tunneling of magneti- 
zation. In the following, we review some of our theoretical work on quantum spin dynamics 
in molecular magnets. 

Ferromagnetic molecular magnets such as Mni2 and Fes show incoherent tunneling of the 
magnetization |3|-|| and allow one to study the interplay of thermally activated processes 
and quantum tunneling. The spin tunneling leads to two effects. Firstly, the magnetization 

1 



relaxation is accelerated whenever spin states of opposite direction become degenerate due to 
the variation of the external longitudinal magnetic field |7HTT|. Secondly, the spin acquires a 
Berry phase during the tunneling process, which leads to oscillations of the tunnel splitting 
as a function of the external transverse magnetic field [p^-^5| . 

Due to the strong quantum spin dynamics induced by antiferromagnetic exchange in- 



teraction |116Hl9|], antiferromagnetic molecular magnets such as ferric wheels belong to the 
most promising candidates for the observation of coherent quantum tunneling on the meso- 
scopic scale pO|-p^]. In contrast to incoherent tunneling, in quantum coherent tunneling 



spins tunnel back and forth between energetically degenerate configurations at a tunneling 
rate which is large compared to the decoherence rate. The detection of coherent quantum 
tunneling is more challenging in antiferromagnetic molecular magnets than in ferromagnetic 
systems, but is feasible with present day experimental techniques. 

Understanding the properties of molecular magnets is only a first step on the roadway 
to technological applications. A possible next step will be the preparation and control of 
a well defined single-spin quantum state of a molecular cluster. Although challenging, this 
task appears feasible with present day experiments and would allow one to carry out quan- 
tum computing with molecular magnets M. The idea is to use the Grover quantum search 



algorithm |24J] to read-in and decode information stored in the phases of a single-spin state. 



Spin tunneling in Mn 12 -acetate 

The magnetization relaxation of crystals and powders made of molecular magnets Mn 12 has 



attracted much recent interest since several experiments [p5|-|29|| have indicated unusually 



long relaxation times as well as increased relaxation rates [fF|,p|,|30| whenever two spin states 
become degenerate in response to a varying longitudinal magnetic field H z . According 
to earlier suggestions pl|j27[1 this phenomenon has been interpreted as a manifestation of 



incoherent macroscopic quantum tunneling (MQT) of the spin. 

As long as the external magnetic field H z is much smaller than the internal exchange 
interactions between the Mn ions of the Mni2 cluster, the Mni2 cluster behaves like a large 



single spin S of length s = 10. For temperatures T > 1 K its spin dynamics can be described 
by a spin Hamiltonian of form 7i = 7i a + Hz + H sp + Ht including the coupling between 
this large spin and the phonons in the crystal [p|- pJ]j3^ -p7[ . In particular, 



H a = -AS* - BS 4 Z (1) 

represents the magnetic anisotropy where A ^> B > 0. The Zeeman term through which 
the external magnetic field H z couples to the spin S is given by Tiz = g^BH z S z , while the 
tunneling between S^-states is governed by 

n T = ~B± {S% + St) + gfi B H x S x , (2) 

where H x = |H| sin# (<C H z ) is the transverse field, with 9 being the misalignment angle. 
The values of the anisotropy constants A, B, B^ have been determined by ESR experiments 
38| , |39[| . Finally, the most general spin-phonon coupling reads 



= 9i(e xx - e yy ) <S> {Si - Sy) + -g 2 e X y ® {S x , S y } (3) 
+ \j9?X^xz ® {S x , S z } + e yz <g> {S y , S z }) + ~gi(u; xz <g> {S x , S z } + u yz <g> {S y , S z }) , 

where ft are the spin-phonon coupling constants, and e a p (u} a p) is the (anti-)symmetric part 
of the strain tensor. From the comparison between experimental data [|7||| and calculation 
it turns out that the constants ft « A Vz P-p[|- 

We denote by \m), —s<m<s, the eigenstate of the unperturbed Hamiltonian TC^ + TCz 
with eigenvalue e m = —Am 2 — Bm A +gfiBH z m. If the external magnetic field H z is increased 
one gets doubly degenerate spin states whenever a level m coincides with a level m! on the 
opposite side of the potential barrier. The resonance condition for double degeneracy, i.e. 
£ m = £m'> leads to the resonance field 

AW = JLu +B t m 2 + m n\\ (4) 

g^B 

As usual, we refer to n = m + m! = even (odd) as even (odd) resonances. 

The relaxation of the magnetization is described in terms of a generalized master equation 
for the reduced density matrix p(t) which includes off-diagonal terms due to resonances ||[10 



We use the notation p mm i = (m \p\ m'), p m = {m \p\ m). In the stationary limit p Tl 
complete master equation 



Pm = -W m p m + ^ W mnPn + ^m (pm' ~ Pm) (5) 

n^m,m' 



can be derived, where 



pm' = F 2 W m + W m , 

m mm 1 . ,„ , , 2 \ / 



is the incoherent tunneling rate from m to mf in the presence of phonon-damping |lOl . The 



spin-phonon rates W m ±i jm and W m ±2,m are evaluated by means of Fermi's golden rule |T(| . 



In Refs. |p|,pro[] the master equation is solved exactly to find the largest relaxation time. 
The result is plotted in Fig. |I[ The even resonances are induced by the quartic i?4-anisotropy, 
whereas the odd resonances are induced by product-combinations of B^S^- and if^S^-terms. 
In Fig. the peaks of the resonance at H z = (induced only by the 54-term) are displayed 
for four different temperatures. All the peaks are of single Lorentzian shape as a result of 
the 2-state transition rate T™' given in (J^j, which agrees well with the measurements [0. 

It is instructive to determine the dominant transition paths via which the spin can relax. 
For this an approximate analytic expression for the relaxation time can be derived by means 
of conservation laws that resemble Kirchhoff 's rules for electrical circuits §JT3 • 



Incoherent Zener tunneling in Fe 8 

Besides Mn 12 there have been several experiments on the molecular magnet Fe§ that re- 



vealed macroscopic quantum tunneling of the spin 1 40- 42, 12, 13 1. In particular, recent mea- 



surements on Fes [fl^ , p~3|1 lead to the development of the concept of the incoherent Zener 
tunneling [|TJJ]. The resulting Zener tunneling probability P; nc exhibits Berry phase oscilla- 
tions as a function of the external transverse field H r . 



For many physical systems the Landau-Zener model jKJ has become an important tool 
for studying tunneling transitions ]4"l]-fi7[] . It must be noted that all quantum systems to 
which the Zener model [^3| is applicable can be described by pure states and their coherent 



time evolution. Ref. [14j generalizes the Zener theory in the sense that also the incoherent 



evolution of mixed states is taken into account (see also Refs. [|3J and |48|-|52|| for a compar- 
ison). In particular, the theory presented in Ref. |T3J agrees well with recent measurements 



of P mc (H x ) for various temperatures in Fes [|12],p~3 



For the Zener transition usually only the asymptotic limit is of interest. Therefore it 
is required that the range over which e mm ' (t) = e m — e m i is swept is much larger than the 
tunnel splitting E mm / and the decoherence rate hj mm ' (see below and Fig.||). In addition, the 
evolution of the spin system is restricted to times t that are much longer than the decoherence 
time Td = l/7 mm /. In this case, tunneling transitions between pairs of degenerate excited 
states are incoherent. This tunneling is only observable if the temperature T is kept well 
below the activation energy of the potential barrier. Accordingly, one is interested only in 
times t that are larger than the relaxation times of the excited states. Thus, the formalism 
presented in Refs. can be applied. It was shown in Ref. ]14j that the Zener tunneling 

can be described by Eq. (||), where 



r™'(t) 



El 



Ir, 



(7) 



2 e mm> W + ^ 2 7mm' 

is time-dependent, in contrast to Eq. (jg). As usual, the abbreviations 7 mm ' = (W m + W m i)/2 
and W m = W nm are used, where W nm denotes the approximately time-independent 
transition rate from \m) to \n), which can be obtained via Fermi's golden rule PJToH. The 
tunnel splitting P,[lO| is given by 

N-l 
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i=l £m ~ em i + l 
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m,N,m' 



(8) 



V m ^ rrij denote off-diagonal matrix elements of the total Hamiltonian 7i 



tot- 



Since all resonances n lead to similar results, Eq. (|^) is solved only in the unbiased 
case — corresponding to n = (see below) — where the ground states \s), \—s) and the 
excited states |m), |— m), m e [[s] — s + 1, s — 1] of the spin system with spin s are pairwise 
degenerate. In addition, it is assumed that the excited states are already in their stationary 
state, i.e., p m = Vm ^ s, —s. Eq. (|5]) leads then to 



1 - P inc = Ap(t) = exp |- J dt' r tot (t')| , (9) 

where Ap{t) = p s — p_ s , which satisfies the initial condition Ap(t = t ) = 1, and thus 
Pmc{t = ^o) = 0. The total time-dependent relaxation rate is given by r tot = 2 [Tg S + T t h], 
where the thermal rate T t h, which determines the incoherent relaxation via the excited 
states, is evaluated by means of relaxation diagrams PUlOfl. 

Assuming linear time dependence, i.e., £ mm >(t) = in the transition region f43|, and 
with \£mm'\ ^ hlrnm' one obtains from Eq. (|) 

Ap = exp i — s ' s arctan ( — — — t J — [ dt' T th \ 
{ ha s s \nis,-s J J-t J 

™^{-il^-[f T A' (io) 

where to = — t. In the low-temperature limit T — > the excited states are not populated 
anymore and thus r\h, which consists of intermediate rates that are weighted by Boltzmann 
factors b m PJlOfl, vanishes. Consequently, Eq. (|l^) simplifies to 



-8, — S 



The exponent in Eq. flTT| ) differs by a factor of 2 from the Zener exponent fI5 |. This is not 



surprising since r to t is the relaxation rate of Ap, where both p s and p- s are changed in time 
by the same amount, and not an escape rate like in the case of coherent Zener transition, 
where only the population of the inital state is changed in time. Eq. flTTD implies P inc = 1 
for |e: Sj _ s (0)| — > (adiabatic limit) and P inc = for |e gj _ s (0)| — > oo (sudden limit). 

In accordance with earlier work [^2|p!3| , [j0| - |42| . |53| Ref. uses a single-spin Hamiltonian 



7i = Ti & + Ht + Hz + Ti-sp that describes sufficiently well the behavior of the large spin S 
with s = 10 of a Feg cluster. After fitting the parameters the incoherent Zener theory is in 



excellent agreement with experiments [|T^,|T5] for the temperature range 0.05 K< T < 0.7 K 
if the states |±10), |±9), and |±8) are taken into account. In particular, the path leading 
through | ±8) gives a non-negligible contribution for T > 0.6 K. Solving the relaxation 
diagram shown in Fig. f| one obtains from Eq. (|To|) for Fes m the case n = 
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tot 



2 1 r r „ 10 



Ap = exp 



n=9 Wio,' 
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(12) 



-kEI _ n W 10 , n b n 



10,n 



where the approximation 7„_ n ~ P^io.n and S> E n _ n ,^ n _ n is used. P ir 

which is plotted in Fig. |5|, is in good agreement with the measurements fl3| . 



1- Ap, 



Coherent Neel vector tunneling in antiferromagnetic molecular wheels 

Antiferromagnetic molecular clusters belong to the most promising candidates for the ob- 



servation of coherent quantum tunneling on the mesoscopic scale currently available [20 



Several systems in which an even number N of antiferromagnetically coupled ions is ar- 



ranged on a ring have been synthesized to date |54]-[5^]]. These systems are well described 
by the spin Hamiltonian 

N N N 

H = J Si ■ s i+ i + gfi B B -^Si-kzY^ ( 13 ) 

i=l i=l i=l 

where Sj is the spin operator at site i with spin quantum number s, s^+i = Si, J is the 
nearest-neighbor exchange, B the magnetic field, and k z > the single-ion anisotropy di- 
rected along the ring axis. The parameters J and k z have been well established both for 
various ferric wheels @-|56|,^8|-^] with N = Q,8, 10, and, more recently, also for a Cr 



wheel ||57|| . For B = 0, the classical ground-state spin configuration of the wheel shows 
alternating (Neel) order with the spins pointing along ±e 2 . The two states with the Neel 
vector n along ±e 2 [Fig. |6|, labeled | |) and | j), are energetically degenerate and separated 
by an energy barrier of height Nk z s 2 . Because antiferromagnetic exchange induces dynamics 
of Neel ordered spins, the states | |) and | j) are not energy eigenstates. Rather, a molecule 
prepared in spin state | f) would tunnel coherently between | f) and | I) at a rate A/h, 



where A is the tunnel splitting [|16],[I7|. This tunneling of the Neel vector corresponds to a 
simultaneous tunneling of all iV spins within the wheel through a potential barrier governed 
by the easy-axis anisotropy. Within the framework of coherent state spin path integrals, 



an explicit expression for the tunnel splitting A as a function of the magnetic field B has 
been derived ||20|| . A magnetic field applied in the ring plane, B x , gives rise to a Berry phase 
acquired by the spins during tunneling p~5| , |63| . The resulting interference of different tun- 
neling paths leads to a sinusoidal dependence of A on 5^ which allows one to continuously 
tune the tunnel splitting from to a maximum value which is of order of some Kelvin for 
the antiferromagnetic wheels synthesized to date. 

The tunnel splitting A also enters the energy spectrum of the antiferromagnetic wheel 
as level spacing between the ground and first excited state. Thus, A can be experimen- 
tally determined from various quantities such as magnetization, static susceptibility, and 
specific heat. Even more information on the physical properties of antiferromagnetic wheels 
[Eq. (0)] can be obtained from a theoretical and experimental investigation of dynamical 
quantities, such as the correlation functions of the total spin S = J2?=i &i or °f single spins 



within the wheel pl|- [23|| . By symmetry arguments, it follows that the correlation function 
of total spin, (S a (t)S a (0)), which is experimentally accessible via measurement of the alter- 
nating current (AC) susceptibility does not contain a component which oscillates with the 
tunnel frequency A/h [[2"T|-|2"B|] . Hence, neither the tunnel splitting nor the decoherence rate 
of Neel vector tunneling can be obtained by experimental techniques which couple to the 
total spin of the wheel. In contrast, the correlation function of a single spin 

/ p -0A/2 p /3A/2 \ 

<m»mo)> - s> (^W a,/ " + ^miTf'n (14) 

exhibits the time dependence characteristic of coherent tunneling of the quantity Sj j2 with a 
tunneling rate A/h ||21|| . We conclude that local spin probes are required for the observation 
of the Neel vector dynamics. Nuclear spins which couple (predominantly) to a given single 



spin Sj are ideal candidates for such probes [21] and have already been used to study spin 



cross-relaxation between electron and nuclear spins in ferric wheels [|64],|65|. 

For simplicity, we consider a single nuclear spin I, I — 1/2, coupled to one electron 
spin by a hyperfine contact interaction H' = Asi ■ I. According to Eq. (|Hp, the tunneling 
electron spin §i produces a rapidly oscillating hyperfine field Ascos(Atfh) at the site of the 
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nucleus. Signatures of the coherent electron spin tunneling can thus also be found in the 
nuclear susceptibility. For a static magnetic field applied in the plane of the ring, B x , it can 
be shown that the nuclear susceptibility 



tanh ( ) 6(u> - HBJH) (15) 



+ (^Ytanh (^)5(u-A/h) 



\LJ — > — U0\ 



exhibits a satellite resonance at the tunnel splitting A of the electron spin system ||21|| . Here, 
77-Ba; is the Larmor frequency of the nuclear spin and the first term in Eq. flTBj ) corresponds 
to the transition between the Zeeman-split energy levels of /. Because typically As ~ 1 mK 
and A < 2 K in Feio, the spectral weight of the satellite peak is small compared to the one 
of the first term in Eq. (|TJ) unless the magnetic field is tuned such that A is significantly 
reduced compared to its maximum value. The observation of the satellite peak in x" yyi^) is 
challenging, but possible with current experimental techniques [[H]]. The experiment must 
be conducted with single crystals of an antiferromagnetic molecular wheel with sufficiently 
large anisotropy k z > 2J/(Ns) 2 at high, tunable fields (10 T) and low temperatures (2 K). 
Moreover, because the tunnel splitting A(B) depends sensitively on the relative orientation 



of B and the easy axis [p9| , p0|| , careful field sweeps are necessary to ensure that the satellite 
peak in Eq. (|T5| ) has a large spectral weight. 

The need for local spin probes such as NMR or inelastic neutron scattering to detect 
coherent Neel vector tunneling can be traced back to the translation symmetry of the spin 
Hamiltonian H [ p2"| , f23f . If this symmetry is broken, e.g. by doping of the wheel, ESR also 
provides an adequate technique for the detection of coherent Neel vector tunneling. If one of 
the original Fe or Cr ions of the wheel with spin s = 5/2ors = 3/2, respectively, is replaced 
by an ion with different spin s' ^ s, this will in general also result in a different exchange 
constant J' and single ion anisotropy k' z at the dopand site, i.e., 

N-l N N 

H = J ^ s-i ■ Si+i + J'i&i • s 2 + si • s N ) + g/j, B B ■ - (k' z s\ z + k z ^ s\z)- ( 16 ) 

i=2 i=l i=2 
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Although thermodynamic quantities, such as magnetization, of the doped wheel may differ 
significantly from the ones of the undoped wheel, the picture of spin tunneling in antiferro- 
magnetic molecular systems |IB|,|I7| remains valid p2| . However, due to unequal sublattice 
spins, a net total spin remains even in the Neel ordered state of the doped wheel [Fig. [7| 
which allows one to distinguish the configurations sketched in Fig. |^ according to their total 
spin. The dynamics of the total spin S is coupled to the one of the Neel vector, and coherent 
tunneling of the Neel vector results in a coherent oscillation of the total spin such that the 
tunneling dynamics can also be probed by ESR. The AC susceptibility shows a resonance 
peak at the tunnel splitting A, 

x'Ulu ~ A/h) = 7r(^ / i i? ) 2 |( e |^|( 7 )| 2 tanh ( ^] 5(uj - A/H). (17) 



2 

with a transition matrix element between the ground state \g) and first excited state |e), 

(18) 



for g^isBx ^> s\/&Jk z . The matrix element in Eq. ( |f8| ) determines the spectral weight of 
the absorption peak in the ESR spectrum. The analytical dependence has been determined 
within a semiclassical framework and is in good agreement with numerical results obtained 
from exact diagonalization of small systems [Fig. [7j . 

In conclusion, several antiferromagnetic molecular wheels synthesized recently are 
promising candidates for the observation of coherent Neel vector tunneling. Although the 
observation of this phenomenon is experimentally challenging, nuclear magnetic resonance, 
inelastic neutron scattering, and ESR on doped wheels are adequate experimental tech- 
niques. The theory of coherent spin quantum tunneling as presented above applies to zero- 
dimensional systems, such as small ferric wheels. With increasing wheel size, the possibility 
of different magnetic domains arises and new exciting quantum effects in the dynamics of 
domain walls come into play f66H?2f. Molecular wheels will also allow one to study the 



transition from zero- to one- dimensional quantum behavior with increasing system size. 



10 



Quantum computing with molecular magnets 

Shor and Grover demonstrated that a quantum computer can outperform any classical com- 
puter in factoring numbers |73[ and in searching a database by exploiting the parallelism 
of quantum mechanics. Recently, the latter has been successfully implemented [71] using 
Rydberg atoms. In Ref. [[]J an implementation of Grover's algorithm was proposed that 
uses molecular magnets [|^,^, |T3iPl , [75| . It was shown theoretically that molecular magnets 
can be used to build dense and efficient memory devices based on the Grover algorithm. 
In particular, one single crystal can serve as a storage unit of a dynamic random access 
memory device. Fast electron spin resonance pulses can be used to decode and read out 
stored numbers of up to 10 5 , with access times as short as 10 -10 seconds. This proposal 
should be feasible using the molecular magnets Fe 8 and Mn 12 . 

Suppose we want to find a phone number in a phone book consisting of N = 2™ entries. 
Usually it takes N/2 queries on average to be successful. Even if the N entries were encoded 
binary, a classical computer would need approximately log 2 N queries to find the desired 
phone number |24j|. But the computational parallelism provided by the superposition and 



interference of quantum states enables the Grover algorithm to reduce the search to one single 
query [J24|. This query can be implemented in terms of a unitary transformation applied to 



the single spin of a molecular magnet. Such molecular magnets, forming identical and largely 
independent units, are embedded in a single crystal so that the ensemble nature of such a 
crystal provides a natural amplification of the magnetic moment of a single spin. However, 
for the Grover algorithm to succeed, it is necessary to find ways to generate arbitrary 
superpositions of spin eigenstates. For spins larger than 1/2 this turns out to be a highly 
non-trivial task as spin excitations induced by magnetic dipole transitions in conventional 
electron spin resonance (ESR) can change the magnetic quantum number m by only ±1. 
To circumvent such physical limitations it was proposed to use multifrequency coherent 
magnetic radiation that allows the controlled generation of arbitrary spin superpositions. 
In particular, it was shown that by means of advanced ESR techniques it is possible to 
coherently populate and manipulate many spin states simultaneously by applying one single 

11 



pulse of a magnetic a.c. field containing an appropriate number of matched frequencies. This 
a.c. field creates a nonlinear response of the magnet via multiphoton absorption processes 
involving particular sequences of a and tt photons which allows the encoding and, similarly, 
the decoding of states. Finally, the subsequent read-out of the decoded quantum state can 
be achieved by means of pulsed ESR techniques. These exploit the non-equidistance of 
energy levels which is typical of molecular magnets. 

Molecular magnets have the important advantage that they can be grown naturally as 
single crystals of up to 10-100 /im length containing about 10 12 to 10 15 (largely) independent 
units so that only minimal sample preparation is required. The molecular magnets are 
described by a single-spin Hamiltonian of the form TC Bp i n = 7Y a + V + H sp + 7Y T P, |l0| , [l4 



where H a = —AS 2 — BS^ represents the magnetic anisotropy (A B > 0). The Zeeman 
term V = g/i^H-S describes the coupling between the external magnetic field H and the spin 
S of length s. The calculational states are given by the 2s + 1 eigenstates of H a + gfisH z S z 
with eigenenergies e m = —Am 2 — Bm 4l + gnBH z m, —s<m<s. The corresponding classical 
anisotropy potential energy E(8) = —As cos 2 9 — Bs cos 4 9 + gfisH z s cos 9, is obtained by 
the substitution S z = s cos 9, where 9 is the polar spherical angle. We have introduced the 
notation m, m' = m — m'. By applying a bias field H z such that gfisH z > E mm >, tunneling 
can be completely suppressed and thus TLt can be neglected |9| JlO|JT^1 . For temperatures 



of below 1 K transitions due to spin-phonon interactions (7"4p) can also be neglected. In 
this regime, the level lifetime in Fe§ and Mn 12 is estimated to be about = 10~ 7 s, limited 
mainly by hyperfine and/or dipolar interactions [|TJ. 

Since the Grover algorithm requires that all the transition probabilites are almost the 
same, Ref. |l| proposes that all the transition amplitudes between the states |s) and |m), 
m = 1,2,..., s— 1, are of the same order in perturbation V. This allows us to use perturbation 
theory. A different approach uses the magnetic field amplitudes to adjust the appropriate 
transition amplitudes ||76|| . Both methods work only if the energy levels are not equidistant, 



which is typically the case in molecular magnets owing to anisotropies. In general, if we 
choose to work with the states m = mo, mo + 1, s — 1, where mo = 1, 2, s — 1, we have 
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to go up to nth order in perturbation, where n = s — mo is the number of computational 
states used for the Grover search algorithm (see below), to obtain the first non- vanishing 
contribution. Fig. [5] shows the transitions for s = 10 and mo = 5. The rath-order transitions 
correspond to the nonlinear response of the spin system to strong magnetic fields. Thus, a 
coherent magnetic pulse of duration T is needed with a discrete frequency spectrum {u m }, 
say, for Mni2 between 20 and 300 GHz and a single low-frequency around 100 MHz. The 
low-frequency field H z (t) = H (t) cos(ujot)e z , applied along the easy-axis, couples to the spin 
of the molecular magnet through the Hamiltonian 

Mow = gfiBH (t) cos(u t)S z , (19) 

where hu <C e mo —e mQ+ i and e z is the unit vector pointing along the z axis. The 7r photons 
of V\ ow supply the necessary energy for the resonance condition (see below) . They give rise 
to virtual transitions with Am = 0, that is, they do not transfer any angular momentum, 
see Fig. |]. The perturbation Hamiltonian for the high-frequency transitions from \s) to 
virtual states that are just below |m), m = m , . . . ,s — 1, given by the transverse fields 
H I(*) = Emimo H m (t)[cos(u m t + $ m )e x - sm(u m t + $ m )e y ], reads 

s-1 

Vhigh(*) = 9^BH m {t) [cos(^ m t + ® m )S x - sm(u m t + ® m )S y ] 

m=mo 
m=mo 

with phases $ m (see below), where we have introduced the unit vectors e x and e y pointing 
along the x and y axis, respectively. These transverse fields rotate clockwise and thus 
produce left circularly polarized o~ photons which induce only transitions in the left well (see 
Fig. |]). In general, absorption (emission) of a" photons gives rise to Am = —1 (Am = +1) 
transitions, and vice versa in the case of cr + photons. Anti-clockwise rotating magnetic 
fields of the form H|(t) = J2m=rn H m (t)[cos(u m t + <3> m )e x + sin(u; m t + ® m )e y ] can be used 
to induce spin transitions only in the right well (see Fig. |[)). In this way, both wells can be 
accessed independently. 
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Next we calculate the quantum amplitudes for the transitions induced by the magnetic 
a.c. fields (see Fig. |9[) by evaluating the S-matrix perturbatively. The jth-order term of 
the perturbation series of the S'-matrix in powers of the total perturbation Hamiltonian 
V(t) = VJow(X) + ^high(^) is expressed by 

j J-i 



(i \ J poo roc 

tz n / / dt 3 e(t k -t k+1 ) 



x[/(oo, ti)V(ti)tf(ti, t 2 ) V(t 2 ) • • • VfcMtj, -oo), (21) 

which corresponds to the sum over all Feynman diagrams T of order j, and where U(t, t ) = 
e -t(?U+gnBSH x )(t-to)/n j g f ree propagator, 0(t) is the Heavyside function. The total S'- 
matrix is then given by S = Y^jLo^ • ^he high-frequency virtual transition changing m 
from s to s — 1 is induced by the frequency lo s -\ = uj s -i iS — (n — l)a>o. The other high 
frequencies u m , m = mo, . . . , s — 2, of the high-frequency fields if m mismatch the level 
separations by uq, that is, hu> m = e m — e m+ i + Jujq, see Fig. ||. As the levels are not 
equidistant, it is possible to choose the low and high frequencies in such a way that S^.s = 
for j < n, in which case the resonance condition is not satisfied, that is, energy is not 
conserved. In addition, the higher-order amplitudes \Sm, s \ are negligible compared to \Sm, a \ 
for j > n. Using rectangular pulse shapes, H k (t) = H k , if — T/2 < t < T/2, and otherwise, 
for k = and k > m , one obtains (m > mo) 

<?<»)- X^O 2n ( 9 » B Y U S k = m H k e-^Hr mo p m ,s^) 
m ' s V m i \2h) (-l^Ir.^H- 1 

x5 {T) (uj mjS - ^uj k - (m -moVo), (22) 



where f2 m = (m — mo)!, is the symmetry factor of the Feynman diagrams (see Fig. [|), 
gjr = m - m - r s (.F), Pm^J 7 ) = U. k=m ( k \ s *\ k) rk(:F \ r^J 7 ) = 0, 1, 2, . . . , < m - m is the 
number of 7r transitions directly above or below the state depending on the particular 
Feynman diagram J 7 , and 5^ t \lj) = ^ f^!^ e^dt = sin(u;T/2)/7ru; is the delta-function 
of width 1/T, ensuring overall energy conservation for ujT ^> 1. The duration T of the 
magnetic pulses must be shorter than the lifetimes of the states |m) (see Fig. [|). 
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In general the magnetic field amplitudes Hk must be chosen in such a way that perturba- 
tion theory is still valid and the transition probabilities are almost equal, which is required 
by the Grover algorithm. In Ref. [jXJ] the amplitudes Hk do not differ too much from each 
other due to the partial destructive interference of the different transition diagrams shown 
in Fig. Ref. |7B[] shows that the transition probabilities can be increased by increasing 



both the magnetic field amplitudes and the detuning energies under the condition that the 
magnetic field amplitudes remain smaller than the detuning energies. In this way, both 
high multiphoton Rabi oscillation frequencies and small quantum computation times can be 
attained. This makes both methods [|IT,|76f very robust against decoherence sources. 

In order to perform the Grover algorithm, one needs the relative phases (p m between 
the transition amplitudes Smjs, which is determined by $ m = Y^k=s-i ® k + fmi where $ m 
are the relative phases between the magnetic fields H m (t). In this way, it is possible to 
read-in and decode the desired phases $ m for each state \m). The read-out is performed by 
standard spectroscopy with pulsed ESR, where the circularly polarized radiation can now be 
incoherent because only the absorption intensity of only one pulse is needed. Fig. |IU| shows 
an example where the state |7) is populated within the computational basis (excluding the 
ground state). Then, we would observe only transitions from |7) to |8) at the frequency 
uj = ugj and transitions from |7) to |6) at u — uqj, which uniquely identifies the populated 
level, because the levels are not equidistant. This spectrum identifies all the populated 
states unambiguously. Alternatively, one could measure the magnetization of the sample 
with a high precision (see e.g. Ref. |77l). We emphasize that the entire Grover algorithm 



(read-in, decoding, read-out) requires three subsequent pulses each of duration T with Td > 
T > Uq 1 > u;" 1 > UJ m'mH- This gives a 'clock-speed' of about 10 GHz for Mni 2 , that is, the 
entire process of read-in, decoding, and read-out can be performed within about 10~ 10 s. 

The proposal for implementing Grover's algorithm works not only for molecular magnets 
but for any electron or nuclear spin system with non-equidistant energy levels, as is shown 
in Ref. J7BJ for nuclear spins in GaAs semiconductors. Instead of storing information in the 
phases of the eigenstates \m) Q, in Ref. [76| the eigenenergies of \m) in the generalized 
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rotating frame are used for encoding information. The decoding is performed by bringing 
the delocalized state (1/y/n) \m) into resonance with \m) in the generalized rotating 
frame. Although such spin systems cannot be scaled to arbitrarily large spin s — the larger 
a spin becomes, the faster it decoheres and the more classical its behavior will be — we can 
use such spin systems of given s to great advantage in building dense and highly efficient 
memory devices. 

For a first test of the nonlinear response, one can irradiate the molecular magnet with 
an a.c. field of frequency u s -2,s/2, which gives rise to a two-photon absorption and thus to 
a Rabi oscillation between the states \s) and \s — 2). For stronger magnetic fields it is in 
principle possible to generate superpositions of Rabi oscillations between the states \s) and 



\s — 1), \s) and \s — 2), \s) and \s — 3), and so on (see also Ref. [7C 
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FIGURES 
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FIG. 1. Full line: semilogarithmic plot of calculated relaxation time r as function of magnetic 
field H z at T = 1.9 K. Dots and error bars: data taken from Ref. [Q. 
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FIG. 2. Full lines: semilogarithmic plots of calculated relaxation rate V = 1/r as function of 
H z for the first resonance peak at (a) T = 2.5 K, (b) T = 2.6 K, (c) T = 2.7 K, and (d) T = 2.8 
K. All peaks are of single Lorentzian shape. Dots: data taken from Ref. 0. 




FIG. 3. Energy level crossing diagram for incoherent Zener transitions. Dotted lines: transi- 
tions due to interaction with environment, leading to a linewidth j mm '. 
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FIG. 4. Unbiased (n = 0) relaxation diagram for Fes- Full (dashed) lines: thermal (tunneling) 
transitions. 
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FIG. 5. Zener transition probability P- m c{H x ) for temperatures T =0.7 K, 0.65 K, 0.6 K, 0.55 
K, 0.5 K, 0.45 K, and 0.05 K. The fit agrees well with data (Ref. p| ). Note that P inc is equal to 



2P in Ref. @ 
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FIG. 6. The two degenerate classical ground state spin configurations of an antiferromagnetic 
molecular wheel with easy axis anisotropy. 
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FIG. 7. The doped antiferromagnetic molecular wheel acquires a tracer spin which follows 
the Neel vector dynamics (left panel). Comparison of results obtained for the matrix element 
Kepzlgf)! with a coherent state spin path integral formalism (solid line) and by numerical exact 
diagonalization (symbols) for N = 4, s = 5/2, s' = 2, J' = J, k z = k' z = 0.055J (right panel). 




FIG. 8. Double well potential seen by the spin due to magnetic anisotropies in Mni2- Arrows 
depict transitions between spin eigenstates driven by the external magnetic field H. 
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Energy 




(5) 

FIG. 9. Feynman diagrams T that contribute to Sm's for s = 10 and rriQ = 5 describing 
transitions (of 5th order in V) in the left well of the spin system (see Fig. 0). The solid and dotted 



arrows indicate a and tt transitions governed by Eq. (|20| ) and Eq. (|19D, respectively. We note 



that S, 



0) 

m,s 



for j < n, and Sm] s <C Sm's for j > n. 
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FIG. 10. Read-out of the information by ESR in the linear regime. In this example only the 
state 1 7) of the computational basis is populated. Thus only the transitions marked by solid arrows 
induce a response that can be observed in a ESR spectrum. 
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